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Abstract. - On the example of monolayer manganites we show that the theoretical ideas developed 
long ago along the derivation of the spin t-J model from the Hubbard model are nowadays very 
helpful in strongly correlated oxides with partly filled degenerate orbitals. We analyze a realistic 
orbital t-J model for eg electrons in Lai_a;Sri+2,Mn04 monolayer manganites, and discuss the evo- 
lution of spin and orbital correlations under increasing doping by performing exact diagonalization 
of finite clusters, with electronic kinetic energy determined self-consistently at finite temperature 
by classical Monte Carlo for t2g spins. Several experimental results are reproduced. 



Spin t— J model. — A new perspective to treat strongly correlated electrons w&s opened 
by the derivation of the spin t-J model three decades ago [1] . We argue that this derivation 
opened new routes to understand the complex magnetic and superconducting instabilities in 
the Hubbard model, and proved to be a very useful paradigm in strongly correlated electron 
systems. At that time it was already realized that the Hubbard model, 



where is an electron creation operator in spin state a =t,i at site i and U is local 
Coulomb interaction, is perfectly designed to describe electron localization in a Mott insu- 
lator, when the average filling is n = 1 electron per site, and the Coulomb repulsion is large, 
i.e., [t/U) ^ 1. In this case the low-energy physics is determined by the effective superex- 
change interaction between localized electron states, which is antiferromagnetic (AF) due 
to the Pauli principle. In this limit the energy is determined by the nearest neighbor spin 
correlations and by higher order correction terms in powers of {t/U) [2]. 

However, when a Mott insulator is doped by holes and electron density is reduced to 
n ^ \ — X, the kinetic energy competes with the AF superexchange and may lead to novel 
physical phenomena [3]. Usually one considers then the effective spin t-J model. 
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derived first in Ref. [1] , where Si are spin operators at each site occupied by a single electron, 
and J = At^ /U is the superexchange constant. In contrast to Eq. (1), the kinetic energy term 
oc t contains 'constrained' fermion creation operators, = a\^{\ — rii^) with a = —a, and 
describes processes restricted to the subspace without doubly occupied sites. The lowest 
kinetic energy is obtained when spins are aligned in a ferromagnetic (FM) configuration. 
This competition between the kinetic and magnetic energy is the driving force responsible 
for the evolution of the magnetic phase diagram under increasing hole doping [4] . 

It was also recognized already in 1977 [1] that the faithful low-energy representation of 
the Hubbard model is not the t-J model itself, but rather the t-J model extended by the 
three-site term which contribute to the hole dynamics away from half filling: 

Weff = -tY^ a\^a^„ + 5, • 5, - i J ^ aLa/ca^feff'^fta' 4<T'«j<T' • (3) 

ija {ij) ikj,<7<7' 

In fact, there is no reason to neglect this last term in Eq. (3) as it is of the same order as 

the effective exchange interaction J [1], and plays an important role except when doping 
a; 0. Indeed, the three-site term captures an important part of the dynamics near a Mott 
insulating state and is essential to give the correct increase of the Drude weight with doping 
X [5], while a constant Drude weight follows instead from the t J model. The three-site term 
is also necessary to reproduce the doping dependence of the spectral weight for the lower and 
upper Hubbard bands in the optical spectra, while the t-J model misses most of the spectral 
weight transfer between the Hubbard subbands which occurs due to kinetic processes [6] . For 
the same reason, only the complete form of the effective Hamiltonian Eq. (3) is appropriate 
to analyze yet another type of order in the regime of large U — the superconducting phase. 
When the local pair operators, bij = (aj^ajj^ — a|^j^a||)/\/2, are introduced [7] for bond (ij), 
the effective Hamiltonian (3) describes correctly the dynamics of local pairs in a strongly 
correlated system. The three-site term gives an important contribution which favors the 
superconducting states over the magnetic ones, and could explain the observed increase of 
the critical temperature in the high superconductors [7]. In contrast, when a propagation 
of a single hole in an AF background is analyzed, the t-J model (2) is sufficient and gives 
new features in the electronic structure, such as quasiparticlc band and incoherent part of 
the spectral density, which result from the coupling of a moving hole to spin excitations [8] , 
while the three-site term leads only to quantitative corrections [9] which become important 
only at large J ^ t [10]. Therefore, the t-J model is also extensively used for interpreting 
the data of photoemission experiments for high Tg superconductors at low hole doping [11]. 
In recent years the t-J model continues to play an important role in the theoretical studies 
of high Tc superconductors, and serves to identify novel complex types of order [12]. 

Inter alia, the idea standing behind the t-J model is extremely helpful to describe the 
doped correlated insulators with orbital degrees of freedom. One decade after the discovery 
of high Tc superconductivity the interest in strongly correlated systems has to a large extent 
moved to the systems with partly filled degenerate orbitals, where electron localization and 
magnetic interactions are rich and lead to complex phase diagrams [13]. In the last decade 
of the last century it was also realized that carriers in doped Mott-Hubbard insulators bind 
to d-d excitations which gives rise to rather general and complex models of t-J variety [14] . 
Both the derivation and theoretical treatment of such modes are a challenge in the modern 
theory of strongly correlated electrons. For instance, internal degrees of freedom of a moving 
carrier can change completely the behavior known from the spin t-J model, where a singlet 
hole is responsible for coherent quasiparticles. In fact, a moving carrier with either orbital 
flavor or large spin (or both) leaves behind a path with its own history. Perhaps the simplest 
example of this behavior is provided by a triplet hole moving in an AF background [15], 
with its completely incoherent motion and suppressed tendency towards FM polarons known 
from the singlet case [4]. 

Yet, there are also cases where the quasiparticle states are obtained by a coupling of the 



p-2 



()rl)ila] corrolatious in mouobu'c^r maiigaiiitcs 



moving singlet hole to orbital excitations. The orbital t-J model [16] was introduced in a 
close analogy to Eq. (2), and describes the motion of a hole in a FM plane of LaMnOa. 
The orbital superexchange J follows again from virtual charge excitations in the regime 
of ?7 3> t, which applies to the realistic parameters of LaMnOa [17]. As in the spin t-J 
model [5], the three-site term plays here an important role near the half filling (for weakly 
doped manganites), and is crucial to reproduce the correct behavior of the Drude weight at 
increasing doping x [18]. Nevertheless, wc neglect the three-site terms as they play almost 
no role at high doping x ~ 0.5 which is of interest in the present study. 

For a FM plane it is allowed to decouple spin and orbital operators following the mean 
field procedure, although in general this would lead to incorrect results. In low-spin states 
the spin-orbital entanglement [19] generates novel ground states which do not follow from the 
classical expectations, such as Goodenough-Kanamori rules [20]. The orbital order (00) 
in an undoped insulator is already a challenge as the orbitals arc directional, the orbital 
(superexchange) interactions depend on the bond direction [14], and their symmetry is thus 
cubic, i.e., much lower than the SU(2) symmetry in the spin case. As a result, the coherent 
quasiparticle states for a single hole moving in a Mott insulator with alternating orbital 
(AO) order have somewhat lower weight than in spin case [16], and depend also on the type 
of the occupied orbitals. 

Only in the last decade it has been recognized that the Jahn- Teller (JT) interactions 
and the superexchange which occurs due to Cg electron excitations in the regime of large 
on-site Coulomb interaction U, support each other [17] and both are necessary to explain 
the magnetic and optical properties of imdoped LaMnOs [21]. In doped manganites, like in 
Lai_a;Sra;Mn03, the average density of Cg electrons n = 1 — x \s tuned by doping, and the 
paradigm of the orbital t-J model may help to understand the actual orbital correlations 
in various doping regimes. In the present paper we follow this idea and analyze the orbital 
correlations in monolayer Lai_j;Sri_|_j,Mn04 manganites using the orbital t-J model. The 
realistic model is richer than that of Ref. [16] , and contains also the JT term [22] , as well as 
the crystal field term which follows from the two-dimensional (2D) geometry and lifts the 
degeneracy of Cg orbitals. The orbital t-J model constructed in this way includes all the 
terms which control the magnetic correlations in doped manganites [23] . 

The monolayer manganites are the subject of intense recent experimental research. The 
undoped monolayer LaSrMn04 compound has the same magnetic structure as K2NiF4, with 
a 2D G-type AF (G-AF) order [24]. At a; = 0.5 the AF order is CE-type, as deduced from 
the x-ray and neutron scattering data for the structural and magnetic correlations [25], as 
well as from the neutron measurements of magnetic excitations [26] . One expects that the 
orbital correlations which follow from strong correlations and the JT interactions play an 
important role in stabilizing this phase [27]. Furthermore, unlike in Lai_,j.Sr,j.Mn03, no FM 
metallic phase was found in doped monolayer Lai_a;Sri_|_a,Mn04 compounds, but instead 
short-range magnetic correlations [24, 25] indicate that magnetic interactions are frustrated. 
This behavior is puzzling and will be addressed in the present study. 

Orbital t—J model. — We consider the effective orbital t-J model, 



introduced here by generalizing the one-dimensional (ID) model of Ref. [22] to doped 
monolayer manganites. The first three terms are the same as the ones in the orbital t- 
J model [16], while the other three serve to describe the realistic situation of eg electrons in 
Lai_a;Sri_|_a;Mn04 manganites. The hopping term Ht for monolayer manganites. 



acts in the restricted space without double occupancies, so the operator creates an 
electron in \a) = \x),\z) state at site i only when this site is initially empty, i.e., c-^ = 



n = Ht+Hj + H,+HjT + Hv+ Hj, , 



(4) 





{i3)\\ab 
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cl^{l — riia) with a standing for the orbital favor opposite to a, similar to spin in Eq. (2). 
We use here the conventional Cg orbital basis: 

|^)^-i=(3z2-r2). (6) 

Here t stands for an effective [ddcr) element that originates from two d p transitions involving 
the intermediate oxygen ion. Following the double exchange mechanism [22], the hopping 
in Eq. (5) is modulated by the quantity Uij which depends on the spin state. Here the core 
spins Si are treated as classical vectors of length S, so their state is represented by two polar 
angles per site, and the scalar product on the bond (ij) is given by 

{S,-S,)=S'{2\u,,\'-l) , (7) 

where the actual dependence on the spins configuration gives 

= COS cos + e'^'I'^-'t'^^ sin sin = e*X'^ cos . (8) 

Therefore, in the classical approximation the hopping (5) is parameterized by the angle 9ij 
between the two involved spins at sites i and j, and by the complex phase Xij- 

The supcrexchange in Lai_j:Sri+2;Mn04 is given by a superposition of several terms 
which originate from charge excitations by either Cg or t2g electrons. The superexchange 
terms which follow from charge excitations by eg electrons are included in Hj. These terms 
favor either FM or AF spin order on a considered bond (ij), depending on the pair of 
occupied eg orbitals at both sites. The form of Hj depends on the actual used orbital basis, 
and it is more convenient to use the directional orbitals and the corresponding pseudospin 
operators, 

= ~HT!tV3T-) , (9) 

depend on the bond direction, with the sign — (+) in Eq. (9) corresponding to a (6) axis. 
The operators are defined by orbital T =1/2 pseudospin operators, 

77 = i< = I {fii. - n«) , Tr = \oi = \ (4c,, + c^c, J , (10) 

with two eigenstates of T/, see Eq. (6). 

The form of Hj used below was derived in Ref. [22] from the complete spin-orbital 
model [17]. It focuses on the orbital dynamics in the presence of spin fluctuations which 
tune orbital superexchange interactions. The spin operators on the bonds (Mn^+-Mn^+ or 

Mn^^-Mn"*"*") arc replaced by their expectation values which involve classical parameters 
given in Eq. (8). In the present case of a monolayer one finds the eg superexchange term [22], 

(y>||a6 

- (1- l^^uf)["ic(l-^j) + (l-«i)"K]} > (11) 

where electron number operator fii^ refers in each case to the directional (3z^ — r^)-like 
orbital \Q along a given bond e.g. — orbital along a axis. The last term in Eq. 
(11) stands for the Mn^+-Mn^+ superexchange and favors configurations of an occupied 
directional orbital next to a hole on the bond, where hi is the electron density operator in 
the restricted space. The superexchange constant J = jU in Eq. (11) is determined by 
the lowest excitation energy C7 to the high-spin 5 = 5/2 state of configuration at the 
involved Mn^+ ion [17]. 

The orbital t-J model [16] contains also a uniform crystal field splitting iJ^ > of eg 
orbitals which removes the orbital degeneracy in monolayer manganites, and is expected 
from the structural data [24]. Therefore, we use 

Hz = ^{fiix - fii;,) , (12) 
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where hi^ and hi^ stand for the electron density operators in the two Cg orbitals. Again, 
as in the spin t-J model, these density operators act in the restricted Hilbert space, with 
double occupancies projected out. If > 0, the \z) orbitals are favored as expected for the 
undoped LaSrMn04 [24], and at a; 0.5 doping in Lao.5Sri.5Mn04 [28]. 
The JT term ifjT includes two contributions: 

Hjr = K T^Tj + ^K ^ THl-nk)Tj. (13) 

{tj}\\ab {{ij)}\\ab 

First term stands for the orbital interaction between two Mn'^^ ions on neighboring sites — 
it plays an important role in the undoped LaSrMn04 and at low hole doping. The second 
one was derived from oxygen distortions by Bala, Horsch and Mack [29] and describes a 
second neighbor interaction along the bonds {{ij)) parallel to either a or b axis; the coupling 
constant fixed here at |«; was estimated from the distance dependence of the JT coupling. 
In addition, we also include the intersite Coulomb interaction between Cg electrons, 

Hv = V h,hj (14) 

{ij)\\ab 

- it plays a role in the doped regime, separating eg electrons from each other. 

Finally, the remaining part of the superexchange follows from charge excitations by t2g 
electrons between two Mn ions. Independently of their valence, the t2g orbitals are filled 
by three electrons which form a core spin S — 3/2 due to Hund's exchange Jh- Therefore, 
this part of the superexchange is AF due to the Pauli principle, in analogy to the AF 
superexchange in the Hubbard model with nondegencratc orbitals. One can verify by using 
the realistic parameters which describe the multiplct spectra of Mn ions that the different 
contributions to the AF superexchange are of similar value [30], so we describe the t2g 
superexchange by the Heisenberg Hamiltonian (J' > 0), 

Hj,=J'Y{Si-Sj-S^). (15) 

(ij) 

For convenience, we consider here core spins Si of unit length. Thereby the actual physical 
values of S = 3/2 spins are compensated by a proper increase of J'. The classical state of 
core spins {Si} is given by two polar angles {19^,0^} per site, which determine the energy 
following Eq. (7), and tune the kinetic energy (5) along each bond {ij). 

The ground state of a monolayer obtained at a given doping x is characterized by electron 
density distribution which we quantify by average electron densities {n^^nz} in orbitals 
a = x,z. The competing tendencies between FM and AF spin correlations at short distances, 
coexisting with different types of 00, have been investigated by intersite spin and orbital 
correlations at distance r: 

= ^ E • ^i+r) , Te{f) = ^Y1 {Ti{0)Ti+m) ■ (16) 

i i 

The orbital operators are defined for a particular orbital basis at site i, 

Ti{9) = cose + sine , (17) 

and the pseudospin operators {Tf,Tf} are given by Eqs. (10). For instance, ^ = corre- 
sponds to TlTl_^„^, 6* = 7r/2 - to TfT^_^„^, and 6* = 27r/3 - to T^T^+p, with ( standing for 
the directional orbital along a axis. The orbital correlations expected in undoped mangan- 
ites are of AO type on two sublattices, which suggests that the orbital correlations Te{r) 
defined as in Eq. (17) are predominantly negative for nearest neighbors. These correlations 
were investigated along the (10) and (11) (and equivalent) directions in 2D clusters. 
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Numerical results. — We employed two different numerical methods to investigate 
finite clusters described by the orbital t~J model (4): (i) exact diagonalization at zero 
temperature (T = 0) for fixed core spin configurations (FM, G-AF, C-AF, and CE phase), 
and (m) Markov chain Monte Carlo (MC) at finite temperature T > for core spins [31], 
combined with exact diagonalization for the eg orbital problem. The ground state at T = 
for representative electron fillings was determined by solving the orbital t-J model for several 
possible types of spin order, using -s/S x -s/S and 4x4 clusters with periodic boundary 
conditions. We solved the orbital problem obtained for the selected core spin configuration 
using the Lanczos algorithm. We then determined the global ground state by comparing 
the energies obtained for different magnetic phases. 

At finite temperatures, we investigated the effective orbital t~J model (4) by making 
use of a combination of Markov chain MC algorithm for the core spins [31] with Lanczos 
diagonalization of the orbital problem. For each classical core spin configuration occurring 
in the MC runs, we defined the actual values of classical variables {uy }, and next solved 
the orbital model. In each case we obtained the free energy for that core spin configuration 
from the few lowest eigenstates, which was next used to decide acceptance in the MC runs. 
More details about the calculation method were given in Ref. [32] . 

We begin with the phase diagram obtained for the undoped LaSrMn04. In absence of 
holes the hopping in Eq. (5) is blocked, and the Coulomb interaction V plays no role as 
it simply adds a constant energy for each bond. Therefore, the ground state for finite J 
is determined by only three parameters: E^/J^ i^/t^ and J' /J. The first two parameters 
determine the orbital state, and the last one tunes the spin interactions. The magnetic order 
depends on all these parameters as the superexchange due to eg electrons is intrinsically 
frustrated and contains both AF and FM terms [17], which become active for particular 
pairs of occupied orbitals on a given bond {ij), see Eq. (11). 

Apart from the unphysical regime of negative J' < which enforces the FM order 
independently of the orbital state, the AO order is necessary to stabilize it for J' ~ and 
\Ez\ < 0.5t. This agrees with the Goodenough-Kanamori rules [20], and with the interplay 
between the magnetic and eg orbital correlations studied before in the ID model [22]. The 
FM phase competes with the G-AF order stable for sufficiently large core spin superexchange 
J' [Fig. 1(a)]. However, the transition line depends also on the two remaining parameters: 
E2 and K. The dependence on Ez is stronger — for E^ < —0.5t the AF order is stable already 
at J' = 0, while for E^ > 0-5t the monolayer would be FM, if the AF superexchange J' 




Fig. 1: (Color online) Magnetic and orbital state for the undoped monolayer (x — 0) obtained with 
an \/8 X \/8 cluster: (a) phase diagram with the stability regions of FM, G-AF, and C-AF phases in 
the {Ez, J') plane at T = 0; (b) electron density n^ in out-of-plane \z) orbitals for increasing value 
of the crystal field splitting E^, as obtained for the FM and AF ground states (solid and dashed 
line), as well as the MC data for pt = 100 and two values of J'. Parameters: J = 0.125t, k = 0.2t. 



p-6 



Orbital correlations in monolayer manganites 




between core spins would not help to stabilize the AF order. This behavior is a manifestation 
of the large difference of 9:1 between the AF eg superexchange interaction [33], depending 
on whether Cg electrons are in |a;) or in \z) orbitals, respectively. Altogether, the transition 
line between the FM and G-AF phase in only weakly dependent on Ez in the regime of large 
\Ez\ > t; once all electrons are found in either \x) or \z) ferro orbital (FO) phase, further 
increase of \Ez\ does not change the state. Remarkably, the range of C-AF phase is very 
narrow and vanishes for E^ > 0.3t, indicating that large overlap between \x) orbitals plays 
a role in stabilizing this type of order, which could also be a finite size effect. 

At the expected value of core spin superexchange J' = 0.025t in LaSrMn04, the G-AF 
state is found independently of E^ [Fig. 1(b)]. However, the experimental data suggest that 
\z) orbitals are occupied in LaSrMn04 [24], and we find that this situation is well described 
by Ez ~ t. When temperature increases to f3t = 100, the orbital state almost does not 
change and indeed the MC data for the density Uz follow the line obtained at T = for 
G-AF phase, see Fig. 1(b). Further increase of temperature promotes excitations to ja;) 
orbitals [24], and this experimental finding is also reproduced by the present model [34]. It 
is interesting to remark that the G-AF order is stabilized in this parameter regime by core 
spin AF superexchange J', while the eg part of the superexchange is FM, as the excitations 
from occupied \z) to unoccupied |a;) orbitals dominate, with the hopping element larger by 
•\/3 than the one between two \z) orbitals, which would give instead the AF coupling. 

At J' = the FM state is stable for Ez > 0, but particularly when the electron densities 
in both eg orbitals are similar for Ez ^ 0, and the 00 may form with occupied orbitals 

\9),eA = cos \z), + sin (p) \x),, \9)j^B = cos {}9) |z), - sin (x), , (18) 

on two sublattices A and B. In this case the FM state is robust and rather large J' ~ 0.025i 
is necessary to switch the magnetic correlations to the G-AF state [Fig. 1 (a)] . This is also 
shown by the MC data for the electron density in \z) orbitals Hz, which remains practically 
unchanged from the ground state value still at temperature pt = 100 [Fig. 1(b)]. The 
density Uz increases slower with Ez than in the AF state, as the AO state (18) can form 
only when the electron density in \z) orbitals is similar to that in \x) orbitals. Note also 
that a rapid drop of Uz with decreasing Ez near Ez — ~0.5t in the MC data for J' — 
accompanies the phase transition to the G-AF phase, see phase boundary in Fig. 1(a). 

We emphasize that, similar to LaMnOa [17], the AO state with FM spin order is stabi- 
lized not only by the orbital superexchange terms in Eq. (11), but also by the orbital JT 
interactions (13). This results in the separation of the energy scales at finite k — 0.2t shown 
in Fig. 2 — the magnetic order is lost first in the range of temperature /3t ^ 50, while the 
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Fig. 3: (Color online) Percentage of electrons found in \z) orbitals n^/n as a function of doping x 
obtained by MC on a x a/8 cluster for three parameter sets: (i) — 0.2t, k — 0.2t, J' = 0.025t, 
(ii) variable = ^{1-x), k = 0.2t, J' = 0.05i, and {in) variable E^ = ~ x), k = 0, f = 0.025t 
(the last set of data has also been shown in Ref. [32]). Parameters: J = 0.125t, V = t. 



AO order is more robust and starts to weaken only above (31 10. While at k = the 
orbital correlations decrease somewhat faster with increasing temperature [34] , they are still 
quite pronounced when the magnetic order is lost. This shows that FM correlations would 
be typically lost at lower temperature T ~ 150 K than the orbital correlations, to a large 
extent due to conflicting trends in magnetic interactions which compete with each other. 

In order to understand the evolution of magnetic interactions with increasing hole doping, 
it is instructive to consider first the orbital occupation. The geometry of the monolayer 
compounds suggests positive crystal field splitting E'^ > of eg orbitals, so the \z) orbitals 
which stick out of the ab plane are favored. We have found that for reasonable positive values 
of Ez the majority of electrons is still in \z) orbitals at low doping, see Fig. 3. However, for 
X > 0.2 both Cg orbitals are either nearly equally populated, or the \x) orbitals dominate, 
i.e., rix > Uz- This behavior confirms the trend found before in the 2D model at orbital 
degeneracy [35]. Here we show that the FO order with \x) orbitals occupied cannot form as 
long as crystal field splitting is not lower than ~ 0. This helps to suppress the kinetic 
energy in ah plane, so the double exchange interaction is weak at finite doping. Actually, for 
this reason the FM state cannot form in the range of intermediate doping 0.2 < a; < 0.4 when 
Ez > Q and the core spin superexchange is in the expected range, J' ~ 0.05t. Altogether, we 
argue that the orbital liquid state which forms only in cubic geometry [36] plays an essential 
role in stabilizing the FM state. 

We believe that the crystal field splitting gradually decreases with increasing eg hole 
density, as the distortions of MnOg octahedra are then expected to decrease and could vanish 
in the limit of a; = 1. Qualitatively this situation is modeled by variable Ez = ^(1 — x)t 
in Fig. 3. While about 70% of electrons fill then \z) orbitals at x = 0, the filling of these 
orbitals at a; = 0.5 is still about 35% (for Ez = \t). In contrast, when Ez = (last set 
of data in Fig. 3), the filling of \z) orbitals drops below 20%, being too low to form the 
observed 00 state [28]. Below we analyze the infiuence of the electron density distribution 
and the microscopic parameters on the stability of the CE-AF phase at a; = 0.5. 

The CE phase was observed in a monolayer Lao.5Sri.5Mn04 compound [25], correspond- 
ing to half doping {x — 0.5). This type of order is quite exotic and there were several 
attempts in the literature to identify the leading parameters which stabilize this complex 
type of coexisting charge, orbital, and magnetic order. We have performed extensive ground 
state {T — 0) calculations using i/S x V8 clusters with periodic boundary conditions, which 
are large enough to capture the essential interactions deciding about the nature of the ground 
state [29], as we also concluded recently [32]. 
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Fig. 4: (Color online) Evolution of the magnetic ground state for x = 0.5 doping as obtained 
with an \/8 x \/8 cluster for increasing crystal field splitting of eg orbitals: (a) energies for 
various possible magnetic phases (G-AF, FM, CE-AF, C-AF) for k = 0.2t and J' = 0.025t; (b) 
phase diagram in {Ez, J') plane for k — 0.2t (solid lines) and k, — OAt (dashed lines). Parameters: 
J = 0.125t, V = t. 



In the absence of JT interactions (at k — 0), the CE phase arises as a compromise 
between competing AF and FM interactions, roughly speaking in a range of < J' < O.lt 
and V < l.bt [32]. In contrast to the undoped case, the magnetic interactions at half doping 
are in general anisotropic and the orbital state can be then tuned to support FM spin 
order along a certain direction by double exchange mechanism, while the superexchange is 
predominantly AF in the other direction. In this way the C-AF (studied before in the ladder 
geometry [37]) and CE-AF phase may arise and strongly compete with each other due to 
rather similar kinetic and interaction energies [32], as long as k = 0. The situation changes 
however when the JT interactions are considered, as finite k — 0.2t used for the data of Fig. 
4(a) favors also alternation of two quasi-planar orbitals, 



\+heA^^{\zh + \x),) I-) 



(19) 



on the charge majority sites in the regime of < -B^ < 0.5t, where electron densities in 
\x) and \z) orbital are similar, rather than directional — r^/3y^ — orbitals. For 
these parameters the C-AF phase is destabilized. Thus, for a moderate value of core spin 
superexchange J' = 0.025t, the CE phase competes only with the FM phase, and the 




-0.4 



O.^i * * * * * * * * * 



-V- v -V- -V - V- 



-V - V- -V - V 



-x--x..-:t 



0.1 



0.2 
K/t 



0.3 



0.4 




Fig. 5: (Color online) Magnetic phases for increasing JT interaction k as obtained with an \/8 x \/8 
cluster for x = 0.5: (a) ground state energies of various phases (G-AF, FM, CE-AF, G-AF) for 
J' = 0; (b) phase diagram in (k, J') plane. Parameters: J = 0.125t, Ez = 0, and V = t. 
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Fig. 6: (Color online) Orbital correlations along the diagonal (1, 1) and (1, 1) directions in the CE 
phase {x = 0.5) for J — 0.125t, and for: V = t, n — 0.2t (solid lines), and V = k — (dashed lines). 

latter one takes over for > 0.15i, i.e., when the population of both eg orbitals is almost 
equal and the double exchange energy dominates. This is also seen in the phase diagram 
of Fig. 4(b) — increasing pushes the CE phase to higher values of J' in a range of 
— 0.25i < Ez < 0.5t. Also the transition between the CE and G-AF occurs at somewhat 
higher values of J' when Ez increases in a range of —t<Ez< 0.25t. Quite remarkably, 
when Ez increases beyond these limits, an opposite trend is found, as then only \z) orbitals 
are occupied and the 00 which accompanies the CE phase (19) cannot develop. Therefore, 
the CE phase collapses, and may exist only in a very narrow window of J' values in between 
the FM and G-AF phase. 

Surprisingly, the dependence of all magnetic energies on the value of the JT interaction 
K is quite weak, see Fig. 5. This interaction stabilizes the 00 both in the FM and in 
CE phase, so the energies of these phases decrease with increasing k [Fig. 5(a)]. However, 
the actual boarder lines between these two phases and between the CE and G-AF phase, 
respectively, are determined by the value of J' [Fig. 5(b)], and depend also on the electron 
density distribution over \x) and \z) orbitals, as discussed above. 

We also note that a finite value oiV ^ t \s important and helps to stabilize the checker- 
board charge order coexisting with the 00 in the CE phase. As expected in this state, 
the |-f) and |— ) occupied orbitals (19) alternate on the sites with increased electron density 
when V = t and k = 0.2t (Fig. 6). This 00 is well developed, but still rather far from the 
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Fig. 7: (Color online) Competition between AF and FM spin interaction at a; = 0.75 doping, as 
obtained with a 4 x 4 cluster: (a) phase diagram in (£2, J') plane with FM, C-AF and G-AF 
phases; (b) MC snapshot of the C-AF phase for — t, J' — O.Obt, fit = 100. Other parameters: 
J = 0.125t, K = Q.2t, and V = t, . 
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idealized one, as Cg electrons are partly dclocalizcd. We note that this partial delocalization 
is crucial for the double exchange which favors the FM order along the zigzag chain. 

Although the competition between the AF superexchange and the FM double exchange 
persists also at higher doping x > 0.5, the checkerboard charge ordered state is destabilized, 
and the optimal conditions for the CE phase are not fulfilled. Instead, one finds a robust 
tendency towards the C-AF phase, as we demonstrate by the data obtained with a 4 x 4 
cluster for doping x = 0.75 in Fig. 7(a). As for x = 0.5, the FM phase is stable at 
J' = 0. In this case the density of Cg electrons is so low that the AF superexchange 
between Mn^+ and Mn''"'" ions is rather weak and the C-AF phase, observed experimentally 
in Ndi_a;Sri+,xMn04 [38], is obtained only when J' > 0.02t, i.e., in the expected range 
obtained by analyzing the t2g superexchange. The value of critical J' decreases somewhat 
with increasing E^, as the double exchange weakens when \z) orbitals are occupied for 
positive values of E^. This effect is quite strong for the second phase boundary, and one finds 
that the C-AF phase survives in a considerably broader range of J' when E^ <0. Actually, 
the critical value of J' exchange changes roughly by a factor of three between E^ = —2t and 
Ez = +2t, which originates from the strength of the double exchange proportional to the 
respective hopping elements between two or two \z) orbitals, respectively. 

Discussion and summary. — The present study clarifies that orbital degrees of free- 
dom are of crucial importance for a complete understanding of magnetic correlations in 
monolayer manganites. We treated a realistic orbital t-J model including Coulomb and 
JT interactions and investigated charge, spin and orbital intersite correlations in monolayer 
manganites. The obtained results revealed a close relationship between orbital and magnetic 
order which follows the classical Goodenough-Kanamori rules at a; = [20]. The magnetic 
phases found in different doping regimes, where double exchange also contributes, are in 
accordance with experiments over the entire doping range < .x < 1. 

For the undoped monolayer LaSrMn04 manganite the model predicts G-AF order, when 
the crystal field splitting of eg orbitals E^ t and the core spin superexchange J' ~ 0.03t 
are in the expected range. Due to large crystal field splitting the role played here by the 
intersite orbital interactions oc k is in this case negligible. We emphasize that the correct 
treatment of electron correlation effects within the orbital t-J model [16] is crucial for the 
qualitatively correct description of the magnetic state in LaSrMn04, and quite different 
results are obtained when strong on-site Coulomb repulsion U is neglected [39] . 

Another success of the model is that it predicts the CE phase at half doping with physi- 
cally realistic parameters for layered manganites, i.e., for small t2g superexchange J' ~ 0.03t, 
as deduced [21] from the analysis of exchange constants in LaMnOs, rather than for unreal- 
istically large J' > 0.2t [40]. The orbital interactions promoted by the JT effect play here a 
very important role in stabilizing this phase, as they support precisely the type of the AO 
order on the checkerboard lattice which promotes FM interactions along the zigzag chains by 
double exchange mechanism. The occupied orbitals are closer to the planar — x^ / — 
orbitals [28], rather than to the directional 3x^ — jiy^ — orbitals, which could also be 
explained by a crystal field Ez> ^ [32], as is frequently believed. Also at doping x ~ 0.75, 
the present orbital t-J model predicts the C-AF phase in agreement with experiment [38] . 

Summarizing, the present study shows that the magnetic phase diagrams of monolayer 

manganites demonstrate a competition between different types of order which follows from 
coexisting FM and AF terms in the spin superexchange, and the FM double exchange at 
finite hole doping. At the magnetic transitions the orbital order changes simultaneously, as it 
supports particular types of magnetic order, in agreement with the Goodenough-Kanamori 
rules. We find it quite remarkable that yet another application of the ideas developed about 
thirty years ago along the derivation of the celebrated spin t-J model [1] turned out to 
be so successful in monolayer manganites, and reproduced several generic features observed 
experimentally in this class of strongly correlated electron systems. 
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